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Abstract. We study the spectral decomposition of the Laplacian on a 
family of fractals VS„ that includes the Vicsek set for n — 2, extending 
earlier research on the Sierpinski Gasket. We implement an algorithm 
[24| for spectral decimation of eigenfunctions of the Laplacian, and ex- 
plicitly compute these eigenfunctions and some of their properties. We 
give an algorithm for computing inner products of eigenfunctions. We 
explicitly compute solutions to the heat equation and wave equation for 
Neumann boundary conditions. We study gaps in the ratios of eigenval- 
ues and eigenvalue clusters. We give an explicit formula for the Green's 
function on VS„. Finally, we explain how the spectrum of the Laplacian 
on VS„ converges as n — > oo to the spectrum of the Laplacian on two 
crossed lines (the limit of the sets V5„.) 
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Figure 1. VS2 Figure 2. VS3 



1. Introduction 

Kigami |16] has developed a theory of Laplacians on a class of fractals 
called pcf self-similar fractals. One example, the Sierpinski gasket SQ has 
become the "poster child" for this theory |i22j in the belief that it is the 
simplest nontrivial example. As a result, a lot of very concrete results have 
been obtained for SQ. This paper extends some of these concepts and results 
to a different family of finitely ramified self-similar fractals, the Vicsek sets 
VSn, with n = 2 corresponding to the Vicsek set VS. We also obtain results 
for VS that have no analogs on SQ. 

To review the standard theory, a pcf self-similar fractal V will be a com- 
pact set in the plane, defined as the limit of a sequence of graphs rojTi, . . . 
with vertices Vq C Vi C • • • . The property of self-similarity takes the form of 
a family of mappings from V to itself, {Fi} which are contractive similarities 
and have the property that V = Fo{V) U Fi{V) U • • • U Fk{V). For example, 
the Sierpinski Gasket is defined by three similarities, each of which sends 
the entire set SG to one of its three smaller triangular component copies. 
We refer to the graph at stage m of the approximation as the mth level 
graph approximation. The Vicsek set (specifically the second order Vicsek 
set VS2, but sometimes simply called the Vicsek set) is the fractal defined 
by the similarities Fi : VS2 VS2 

Fi{x)=x/3 ^4(2;) = x/3-h 2/3(0, 1) 

F2{x) = x/3 + 2/3(1, 0) Fo{x) = x/3 + 2/3(1/2, 1/2) 

F3(x) =x/3 + 2/3(1,1) 

The first graph approximation Fq is the complete graph on four vertices 
(that is, the vertices of a unit square and an edge connecting every pair of 
vertices). The next approximation Fi consists of five miniature copies of 
Fq arranged in an X shape with branches of length 2 (hence VS2). Further 
graph approximations likewise consist of five copies of the previous level; 
they display finer levels of branching. Higher order Vicsek sets VSn are 
similar, except that Fi is an X-shaped graph consisting not of five but An — 3 
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copies of To, with arms of length n. Instead of five similarities, we have 4n— 3 
similarities. 

It is intuitive from the picture and also easy to demonstrate that as n ^ 
cxD, VSn approaches the pair of crossed line segments between (0,0) and 
(1,1) and (1,0) and (0,1). (That is, the maximum Euclidean distance of 
any point in V5„ from the crossed lines approaches zero.) This is important 
to note because it suggests a connection between fractal analysis on the 
Vicsek sets and classical analysis on the line; later in this paper we show 
that the spectrum of the (Neumann) Laplacian as defined on the Vicsek sets 
does, in fact, approach the spectrum for the classical Neumann Laplacian 
on the cross. 

On VS n, we can define a standard self-similar probability measure as 
follows: for each graph approximation, let Vm be the probability measure 
which weights each vertex by its degree. Then the standard measure fi on 
VSn is defined by 

/ fdfi = lim / fdum- 

We define the unrenormalized energy of a function on T^ by 
Em{u) = \u{x) - u{y)\'^. 

The renormalization factor for VSn is 2n — 1, so the renormalized graph 
energy on Vm is 

and we can define the fractal energy E{u) — IiiHttt, >.(X) £m{u). We define 

dom £ as the space of continuous functions with finite energy. 

Now we have the tools to define a fractal Laplacian. In domiS, £ extends 
by the polarization formula to a bilinear form £{u, v) which defines an inner 
product in this space. If ^ is the standard measure, we can define the 
Laplacian with a weak formulation: A« = / if / is continuous, u G dom£^, 
and 

£{u,v) = — J fv dfi \/v £ domoif 

where domo<£^ = {v G £ : tJlbdry = 0}. There is also a pointwise formula 
(which is proven to be equivalent in [22]) which, for nonboundary points in 
VSn computes 

Au(x)= lim K(4n-3)'^(2n- l)'"A™u(x), 

m— >oo 

with K a constant, and where is a discrete Laplacian associated with 
the graph L^, defined by 

Amu{x) = y^(ti(y) — u{x)), for x not on the boundary. 

degx ^ 
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The Laplacian satisfies the scahng property 

A{u o Fi) = (4n - 3)(2n - l)(An) o Fi 

and by iteration 

A(u o F^) = {{An - 3)(2n - l))"'(Au) o 

for Fu, = Fu,^o Fui^o ■ ■ ■ o • 

In this paper, we restrict attention to the Laplacian defined with Neumann 
boundary conditions. The Neumann boundary conditions are "natural", in 
the sense that the weak formulation need only be modified to allow all 
V S doniiS, and the pointwise formulation is also valid at boundary points. 
It is also possible to define a normal derivative dnu{qi) at boundary point, 
and the Neumann condition is dnu{qi) = 0. Moreover, there are infinitely 
many points in VSn that have neighborhoods isometric to neighborhoods 
of boundary points; the Neumann boundary conditions treat the boundary 
points no differently from these equivalent points. (Note that this is not true 
on SQ.) These are ample reasons to prefer Neumann to Dirichlet boundary 
conditions. An additional benefit is that the theory is considerably simpler. 

The Laplacian on a fractal such as SG or VSn has a discrete spectrum of 
positive eigenvalues Aq < Ai < A2 < • • • , which can be computed explicitly 
by the method of spectral decimation developed by Fukushima and Shima, 
and applied to the Vicsek set in [23]. Spectral decimation is a method of 
relating eigenfunctions and eigenvalues from one graph approximation to a 
finer one. In Section [2| we describe the method and explicitly compute an 
algorithm for spectral decimation on VS2, which allows us to numerically 
calculate eigenfunctions on the Vicsek set, and observe patterns in the data. 

Let {Aj} denote the spectrum of the Laplacian, and let {uj} denote an 
orthonormal basis of eigenfunctions. Then for any bounded function /, we 
can define the spectral operator /(—A) on L^(V5„) by 

00 

f{-A)u = Y,f{X,){u,u,)u,. 
i=i 

These operators include the fundamental solutions to the heat and wave 
equations, and solutions for other space-time equations. Because of the 
importance of spectral operators to classical analysis, understanding spectral 
operators and the Laplacian on VS is a key goal in the development of 
analysis on fractals. 

In computing a spectral operator, we can group terms in the sum corre- 
sponding to the same eigenvalue, and write 

/(-AMx) = ^/(A) / Px{x,y)u{y)df,{y) 

where, at a given point x 

P\{x,y) = ^Uj{x)uj{y), 
j 



THE LAPLACIAN AND SPECTRAL OPERATORS ON THE VICSEK SET 



5 



{uj} being an orthonormal basis of the A-eigenspace Ex. In Section [s] we 
show how, for certain special points x, we can simplify this sum to a single 
term. Fixing a point x on the boundary, or at the center, and letting Ef^ 
denote the subspace of Ex of functions vanishing at x, we can choose the 
orthonormal basis so that the first element ui is in (E^)-^ and the rest belong 



to E^. Then, 



Px{x,y) = ui{x)ui{y). 



Additionally, in Section [3] we prove a formula for the inner product of two 
eigenfunctions on a graph approximation, and show that it converges in the 
limit to the inner product on the Vicsek set. This ensures that functions 
which are orthogonal on graph approximations remain orthogonal on the 
Vicsek set, and makes it possible to compute P\ when a; is a point on the 
boundary or at the center. Here we follow some of the ideas in [2 . 

In Section |4j we give some numerical data using our MATLAB algorithms 
for the eigenvalues and eigenfunctions of the Laplacian on VS2 and V5„. 
We also give data on the eigenvalue counting function N(x) and the Weyl 
ratio N{x)/x°', for the appropriate power a. 

In Section [5] we give numerical results for the heat kernel, the propagator 
for the wave equation, and the spectral projections onto the 0-series. 

In Section [6j we show that each 0-series eigenfunction is determined by 
its restriction to the diagonal of the Vicsek set. 

In Section [7] we prove, following [5] the existence of a ratio gap in the 
spectrum of the Laplacian. A ratio gap is an interval (a, h) such that the 
ratio of any two eigenvalues must fall outside the interval; this is a measure 
of the sparseness of the spectrum. Related results have been obtained in 

m- 

In Section [8] we show the existence of eigenvalue clusters; that is, arbi- 
trarily many distinct eigenvalues in an arbitrarily small interval. 

In Section [9] we calculate an explicit Green's function for the Laplacian 
on the Vicsek set. 

In Section [Toj we examine the convergence of eigenfunctions and eigen- 
values of the Laplacian on VSn as n ^ 00 and show that they approach the 
corresponding values for the Laplacian on the cross. 



In Section 11 we establish some properties of the Weyl ratio on VtS^ 
that begin to explain the curious apparent convergence to a function that is 
unrelated to the Weyl ratio on the cross. 



For more data and programs, refer to Iwww . math . Cornell . edu/~inliw33 

m)- " 

It is possible to describe VSn as the closure of a countable union of straight 
line segments; start with the two diagonals, and take all images under all 
iterates of {-Fj}. (Some images will be proper subsets of other line segments 
and should be deleted to eliminate redundancy.) We call this the skeleton 
of VSn, SKiySn) = U^ilj, where the line segments Ij intersect only at 
points. Since the skeleton is dense, any continuous function is uniquely 
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determined by its restriction to the skeleton, but the skeleton is not all of 
the Vicsek set, since it has /U-measure zero. 

Each line segment Ij has a simple one-dimensional energy 

Sj{u,v)= I' u'{sj{t))v'{sj{t))dt 
J-bj 

where Sj : [—bj,bj\ — > Ij is the linear parametrization. It is not difficult to 
see that 

oo 

S{u,v) = c^£j{u\i.,v\i.) 
i=i 

for the appropriate constant c. Prom this point of view, the energy form on 
VSn is trivial. Because we combine the trivial energy with the unrelated 
measure fi, we obtain a nontrivial Laplacian. 

On the other hand, there is a natural measure on the skeleton: just take 
the sum of Lebesgue measure on each Ij . By the embedding of the skeleton 
in VSn we may also regard this as a measure u on VSn- Of course it is not 
a finite measure, as the sum of the lengths of the line segments Ij diverges. 
It satisfies the self-similar identity 

2n - 1 ^ ' 
in contrast to the self-similar identity 

4n - 3 ^ ' 

for 11. There is good reason to consider i/ as the universal energy measure 
on VSn- If / G domiS then we may define an associated energy measure Uf 
with £{f,f) = Vf^VSn) and roughly speaking i^f{A) is the contribution to 
£{f,f) coming from the set A, for any simple set A (for example, a finite 
union of cells.) For each Ij consider the function fj defined by fj{sj{t)) = t 
on Ij, which is constant on every other interval that intersects Ij- Then 
fj is harmonic at every point except the endpoints of Ij, and Uf. is exactly 
Lebesgue measure on Ij. So 

oo 

We can also see that / = Z^j^i fj is a finite sum on each Ij and v = Vf, 
although / does not have finite energy. One can also show that Vf <^v for 
every function / G dom£^. This is the "universal" property of v- 

On SQ one can define the Kusuoka measure v = + v^^ where {hi, 
is an orthonormal basis of global harmonic functions (modulo constants) in 
the energy norm, and this serves as a universal energy measure. A sim- 
ilar approach would not work on VSn, since it would produce a measure 
supported on the two diagonals alone. 
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It is possible to define an energy Laplacian on VSn using the energy £ 
and the energy measure v in place of /x, although there are some technical 
problems because ly is not finite. Such a Laplacian would be rather "trivial" , 
since it would amount to the second derivative along each line segment Ij, 
together with matching conditions on first derivatives at points of intersec- 
tion. We will not consider this Laplacian further in this paper. 

We hope that this paper makes a strong case that the Vicsek sets deserve 
to be considered the simplest nontrivial examples of pcf self-similar fractals. 
There are two sides to this statement. The first is that the analysis is 
nontrivial. Indeed, if you just restrict attention to harmonic functions on 
VSn, the theory is basically trivial: these are just linear functions on each 
of the arms of VSn that are constant on all trees that attach to an arm. But 
the graphs we have obtained for eigenfunctions of the Laplacian reveal that 
these are nontrivial functions. 

The other side of our assertion is that VSn is simpler than SQ. The 
expression for the Green's function and the numerical data for solutions of 
the wave equation are good a posteriori evidence for this. We can also point 
to two structural features that can be considered a priori evidence. The 
first is topological: VSn is contractible while SQ has infinite dimensional 
homology. Indeed, the cycles in SG play a role in the description of the 
structure of some of the eigenspaces of the Laplacian (the 5-series in the 
terminology of |22j.) The second relates to symmetry: while SQ only has 
a 6-element symmetry group, VSn has an infinite symmetry group. Indeed 
this group is a semidirect product of one copy of 5*4 and infinitely many 
copies of S3 and 5*2. {Sk denotes the permutation group on k letters.) The 
5*4 symmetries are the permutations of the 4 arms, which fix the center 
point go- For any cell, F^V with center point F^qo, with Wm 7^ 0, there 
will be either S2 or 5*3 symmetries permuting 2 or 3 of the arms of the 
cell, depending on whether the cell Fy^V has 2 or 1 neighboring cells (the 
permutable arms are the ones with no neighbors.) 



2. Spectral Decimation 

The method of spectral decimation was invented by Fukushima and Shima 
[12' for SG to relate eigenfunctions and eigenvalues on the graph approx- 
imations to each other and the eigenfunctions and eigenvalues on SQ. In 
essence, an eigenfunction on Tm with eigenvalue Am can be extended to 
an eigenfunction on Tm+i with eigenvalue Xm+ii where Xm = R{Xm+i) for 
an explicit functions R, except for certain specified forbidden eigenvalues, 
and all eigenfunctions on SG arise as limits of this process starting at some 
level m. This is true regardless of the boundary conditions, but if we spec- 
ify Dirichlet or Neumann boundary conditions we can describe explicitly 
all eigenspaces and their multiplicities. This method was extended to the 
Vicsek sets by Zhou [211 . 
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We describe the procedure briefly here. First, there is a local extension 
algorithm that shows how to uniquely extend an eigenfunction u defined on 
Vm to a function defined on Vm+i such that the A-eigenvalue equations hold 
on all points of Vm+i \ Then there is a rational function R{\) such that 
if u satisfies a Am-eigenvalue equation on T^, then the extended function 
will satisfy the A^+i -eigenvalue equation on Vm+i if Am = R{^m+i) and Am 
is not a forbidden eigenvalue. (Forbidden eigenvalues are singularities of the 
spectral decimation function R. It is "forbidden" to decimate to a forbidden 
eigenvalue. Because forbidden eigenvalues have no predecessor — there is 
no Am-i corresponding to Am — we speak of forbidden eigenvalues being 
"born" at a level of approximation m.) 

We have the following theorem from 



Theorem 2.1. Define 

fn{X) = T„(3A - 1) - 3T„_i(3A - 1) 

5n(A) = C/„-i(3A - 1) - C/„_2(3A - 1) 

hn{X) = Un-l{3X - 1) - 3C/„-2(3A - 1) 

where Tn and Un are the Chebyshev polynomials of the first and second kind. 
Then the spectral decimation function R is 

Ri,XjYi) — Am— 1 — Xm9m{XjYi)hfn(^XYn) . 

Moreover, the forbidden eigenvalues are 4/3 and the zeroes of fn and gn. 

We also have a matrix equation for the eigenfunction extension formula: 
If u|vo is a vector of the values of u on Vq and ti| Vi\vb is defined analogously, 
then 

'^IvAVb = -(^ + XiUy^Juvo- 

where J is the Vq x (Vi \ Vq) adjacency matrix, X is the adjacency matrix 
for Vi \ Vq, with the degrees of each vertex as its diagonal entries, and M is 
a diagonal matrix with Ma = —Xa. Multiplying this matrix by the values 
of u on any /c-cell (with Ai replaced by Afc+i), we similarly get the values of 
u on the {k + l)-cells contained in that A:-cell. 

In the case of VS2, we have R{X) = 36X^ - ASX"^ + 15A. The forbidden 
eigenvalues are 0, 1/2, 4/3, and (7 it -v/T7)/12. There is a 0-eigenvalue born 
at level 0, and a 4/3 eigenvalue born at every level thereafter, and continued 
eigenvalues are formed by successively choosing one of the three inverse 
functions of R (see Figure [s]) , so long as this does not lead to a forbidden 
eigenvalue. Using the labeling system described in Figure |4j the matrix 
which allows us to continue eigenfunctions is given by 



(2.1) -{X + XM)-^J 
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cdcaabdcccdc 
cdcccdbaacdc 
\cdcccddccaba) 
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where 

a = 9-42A + 36A^ c=l, 

6 = 6(1-4A + 3A2), d = 2-3X, 

1 

^ ~ 3(4 - 29A + 60A2 - 36A3) ' 

(Note that only roots of —3/7 = (1 — 2A)/(A) are forbidden eigenvalues, so 
7 is well-defined as long as A is not forbidden.) 

We denote the 4/ 3- series as those eigenvalues continued from a 4/3- 
eigenvalue, and the 0-series as those eigenvalues continued from the 0- 
eigenvalue. To find Am from Am-i we have to invert R; in the case of 
VS2, there are three inverses, shown in Figure [3] Note that for the sequence 
IS^Aj to converge to an eigenvalue A on V5, we need Am to approach zero, 
so we must choose the smallest of the three inverses all but finitely many 
times. 

A proof in [2l] guarantees that spectral decimation produces all possible 
eigenvalues and eigenfunctions (up to linear combination), so this formula 
allows us to explicitly determine the values of eigenfunctions at arbitrarily 
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high graph approximations. We make several observations from numerical 
calculation of the eigenfunctions (see Section [g]) . One is that the restrictions 
of certain eigenfunctions to the diagonal (the segment in between (0, 0) 
and (1, 1)) are periodic with period proportional to 1/m and approximate 
sine functions; this suggests that higher Vicsek sets V5„, as they converge to 
a cross, will have eigenfunctions approaching the sine and cosine functions 



in the classical case. We will prove this fact in Section 10 

Secondly, we observe that for the 0-series eigenfunctions, choosing the 
smallest inverse function of R first means that Ai = so the eigenfunction 
is extended to be constant on Vi. On each of the five 1-cells, we start as be- 
fore, with the eigenfunction having a value of 1 one all boundary points; so 
the eigenfunction is miniaturized into identical copies at each graph approx- 
imation, and the eigenvalue is multiplied by 15. The same thing happens 
for any number of initial choices of the smallest inverse function. 

We next describe the structure of the spectrum of the Neumann Laplacian 
on VSn in complete detail. Let (j)2, ■ ■ ■ , 4>2n-i denote the inverse functions 
of the polynomial R in Theorem 2.1, in increasing order. We note that 
is an increasing function when j is odd and is a decreasing function when j 
is even. We write pn = (4n — 3)(2n — 1) for the Laplacian renormalization 
factor. We write = Aq < Ai < A2 < • • • for the distinct eigenvalues. The 
spectral decimation rules are summarized as follows: 

(i) Each eigenvalue has the form 

Yun^p'^ct)^,^ o (t)^^_^ o . . . o </.^^(0), 

JHSo Pn'^^^^^ ° ° • • • ° 0«,i(4/3). 

where in the first case the eigenvalue is in the 0-series and in the 
second case it is in the 4/3-series and born on level k. 

(ii) All but a finite number of the Wm are equal to 1 . 

(iii) For the 0-series, the first Wj with Wj 7^ 1 must have Wj an odd 
number; for the 4/3-series, wi must be an odd number but wi / 
2n- 1. 

(iv) The multiplicity of each 0-series eigenvalue is 1 , while the multiplicity 
of each 4/3-series eigenvalue born on level Zeis 2{An-2,f + 1. 

Condition (ii) is required in order that the limits in (i) exist. Let uiq denote 
the largest value of m for which Wm 7^ 1 (if this never happens, let ttlq = 0.) 
Then we can rewrite the limits in (i) in terms of a single function ipn defined 
by 

v^„(t) = lim p':4>^r\t) 

m — >oo 

(here <^i"^^ denotes the m-fold composition of 0i). This limit exists because 
the Taylor expansion of R{t) about t = is pnt + OijP'), so the Taylor 
expansion of (/>i(i) about t = is p-H + 0{t^). Then (i) says the eigenvalues 
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Pn>n(</'«;„o O </'«;„o_i O " " " O (0)), 
or pr+'V'n(</'«;^„ o </,^^^_i O . . . <^^,(4/3)). 

Condition (iii) spells out explicitly the rules for avoiding forbidden eigen- 
values. Wc may explain the multiplicities in (iv) as follows. To satisfy the 
4/3-eigenvalue equation on level k we may assign initial values at the points 
in Vk so that the sum of the values on the four boundary points of every 
A;-cell is 0. This gives a space of dimension #{Vk} — #{k-cells} and it is 
easy to see that #(14) = 3(4n - S)*^ + 1 and #{A;-cells} = {An - 3)*. 

Theorem 2.2. Eigenvalues in the 0-series and 4/3-series alternate: Xj is 
0-series for j even and 4/3-series for j odd. More precisely, the spectrum 
consists of an initial segment of length 2n followed by segments of length 
An — 2. In each segment all the 4-/3-series eigenvalues are horn level 
(hence have multiplicity 3) except the last one. 

Proof. Because is increasing, so is ipni so applying ipn does not change the 
order of eigenvalues. Thus the ordering of the eigenvalues can be inferred 
from the ordering and increasing/decreasing nature of the (pj^s. The lowest 
eigenvalues have mo = 0, Aq = 0n(O) = 0, and Ai = il)n{A/'i). After that 
come those with mo = 1, in the order 

pnV'n(03(O)) < PnV'n((^3(4/3)) < p„Vn(<^5(0)) < pn(t>n{4>b{A / ?>)) < ■ ■ 

< PnV'n(02n-l(O)). 

(here we have used the fact that </'2j-i is increasing and max 02^-1 < 
min(^2j+i and (/>2n+i(4/3) is a forbidden eigenvalue.) Together these give 
the first 2n — 1 eigenvalues of the initial segment, with 0-series and 4/3-series 
born at level alternating. The segment is completed by A2n-i = Pn^n(4/3), 
a 4/3-series eigenvalue born on level 1. This is valid because 02n-i(O) < 4/3 
(indeed max^2n-i < 1) and p„Vn(4/3) < /OnV'n(min ^2)- 
Let Si denote the sequence 

<^i(4/3), <^3(0), <^3(4/3), 05(0), 05(4/3), . . . , 02n-i(O), 4/3 

and let Si denote the sequence in reverse order with the first and last 

terms omitted. Then the initial segment of the spectrum has the form 
{0,/J„Vn(Si)} (note that ■i/'n(4/3) = /JnV'n((0i(4/3)). Let S2 denote the 
sequence 



0l(Sl), 02(Si), 03(Si), 04(Si), <?!'5(Sl), . . . , 

02n-2(Si), 02n-l(Sl \ {4/3}), 4/3. 
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Then {0,p^Vn(S2)} is a larger initial segment of the spectrum. Note that 
each of the segments 

has length 4n — 2 and alternates 0-series and 4/3- series, where all except the 
last 4/3-series eigenvalues are born on level 0. 
Inductively, we define S^, to be the sequence 

4>l(X'k-l),4'2{,^k~l),<p3{.^k~l), ■ ■ • ,</'2n-2(Sfc_i),(/>2n-l(Sfc_i \ {4/3}), 4/3. 

Then {0, p^tpnC^k)} is an initial segment of the spectrum, and after {0, pnipni^i)} 
it breaks up into segments of length 4n — 2 with 0-series and 4/3-series al- 
ternating, and all but the last 4/3-series alternating, and all but the last 
4/3-series eigenvalues are born on level 0. □ 

2.1. Scaling Inner Products. In order to find an orthonormal basis for 
eigenspaces, we have to relate the graph inner product (/, g)m to the inner 
product on the next graph approximation, {f,g)m+i- This is necessary 
because we need to compute the inner product exactly, and we would like 
to be able to show that functions orthogonal on one graph approximation 
will remain orthogonal when spectrally decimated at higher levels. We now 
prove, as [18^ does for the Sierpinski Gasket, a multiplicative formula for 
(/, g)m+i in terms of (/, g)m and the current discrete eigenvalue A^. 

Theorem 2.3. // u and v are eigenf unctions born on level ml < m, both 
with the same graph eigenvalues \m~i cLnd Xm, then 

{u,v)m = N{m){u,v)m~i 

where 

20 - 143A^ + 240A^ - 108A^ 
4 - 28A„ + 60A2^ - 36A3, 

The product below converges, and in the limit gives the inner product on 
VS2 for u and v eigenf unctions born on level with the same eigenvalue: 

00 

{u,v) = {u,v)q JI N{m). 

m=l 

Proof. On a graph approximation of the Vicsek Set, we call two points 
neighbors if they are connected by an edge. All points have either three 
or six neighbors. We define junction points to be those with six neighbors, 
and non-junction points to be those with three neighbors. 

For simplicity we take u = v as the general case is essentially the same. 

The graph inner product of two functions on the graph approximation Vm 
is defined as 

\ 

(ti, = I • ^ '^u{Fy,qi)v{Fyjqi), 

\w\=m i 
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where we need to multiply by | so that (1, l)m = 1- This makes the limit fi 
a probability measure. Here each w is a "word," that is, a string of numbers 
corresponding to the five similarities that define V52 • F^, is the composition 
Fi-^ o Fi^ o • • • o Fi^ where ij are the constituents of the word w. 

At each graph approximation, these similarities map two distinct points 
to the junction points, and only one point to the boundary points, so we 
account for double-counting as follows: 

{u,v)m = ^-5'"'(2 u{x)v{x)+ u{x)v{x) 

junction nonjunction 

X 

Fix an (m — l)-cell C and let ui, U2, 1*3, U4 be the values of u on its 
boundary. Then the contribution to ||^i||m-i to C is 

Now we can use spectral decimation (i.e. (2.1) with A = Am) to get the 
values of u on the level m vertices in C. Letting a, b, c, d, and 7 as in (2.2 ) 
we see that the contribution of C to 11 Ul I I ,yy, IS 



4'^ 



{ui + ui + ui + ui){l + 2Y{a' + ¥ + 3c^ + U')) 



+ 87^ Y UiUj{ac + bd + c^ + d^) 

i<j 

Applying this to all (m — l)-cells we obtain 

(2.3) ||n||^ = l\\u\\l_,{l + 272(a2 + 6^ + Sc^ + 3d^)) 

+ 2 • S-'^j^ac + bd + c^ + d^) ^ u{x)u{y), 

x^y 

where the sum is over Vm-i- To deal with the cross-terms in terms we apply 
the Gauss-Green formula, 

Em-i{u) = Y - ^y)? 

x^y 

= -12-5""-^(u,A„„in)„_i 
= 12-5"-iA„_i||n||^^_i. 

Since 

Y - u{y)\'^ = ^(deg j;) u{xf - 2 ^ u{x)u{y), 

xr^y X xr^y 

this implies 



^n(x)n(y)=6-5™-i(l-Am_i) 



l^llm-l- 



xr^y 
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Combining this with (2.3), we see 
1 



(2.4) ||n||^ = 9 (1 + 27^(«^ + + 3c^ + 3d^)) 



2/'„2 I l2 I o„2 I oj2\ 



+ 12 • 7^(1 - A™_i)(ac + 6a! + + d^) 



l7/l|2 



Simphfying using the values for 7, a, 6, c, d, and Ar„_i in terms of A^^, we 
get the normalization formula 

, 1 20 - 143A^ + 240A^ - 108A^ 

5 ' 4 - 28Am + 60A2^ - 36A3^ 

This allows us to compute the norm of a function the Vicsek set at any 
graph approximation, and, in the limit, on the Vicsek set itself: 



00 



|2 



]u\\' = \\u\\q ■ Y[ N{m). 

m=l 

□ 

2.2. Center values. It is also useful to have a formula for the value of an 
eigenfunction at the center qq of VS2- Using (2.1) to continue a function u 
on Vq to Vi, we see that the values uq, uiq, uu and uis are related to the 
values of u on Vq by 

+ uio + nil + 1*15 = 7(3(i + b){ui + U2 + U3 + U4). 

Substituting for d, b, and 7 we get 

4-3Ai , 

ug + uio + Uu + ui5 = ^rr^]X^KL8Af + + U4). 

Continuing this process we get 

f U2 

4 

m=l 

where 

4-3A^ 



u{qo) = [[ N {m) 



N'{m) 



4-21A^ + 18A2^' 

In particular, since any 4/3-series eigenfunction satisfies ui+U2+U3+Ui = 0, 
all 4/3-series eigenfunctions vanish at qq. 

3. Spectral projections at boundary points 

We would like to be able to solve differential equations such as the wave 
equation 

52 



5^ = ^^ 



and the heat equation 

du 
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on the Vicsek set with Neumann boundary conditions and suitable initial 
conditions. These equations are solved in terms of an orthonormal basis of 
eigenfunctions uj: we have 

u{x,t) =J2Sit,j) (^J f{y)uj{y)dfj,{y)^ Uj{x) 

where S depends on the equation and the expression in parentheses is a 
Fourier coefficient. Usually the sum and integral can be interchanged to 
yield 

u{x,t) = J Kt{x,y)f{y)d^i{y), 

where / is defined by the initial conditions and 

Kt{x,y) = ^S{j,t)uj{x)uj{y), 

for instance 

Y,e~'^'Uj{x)uj{y) 

for the heat equation and 

J2 ^^^J^^* Uj{x)uj{y) 

for the wave equation. 

We can get a better understanding of projection kernels Kt{x,y) when 
we restrict one of the variables to specific boundary points. Suppose y = Qi, 
i = 1,2,3,4, and note that E^ = {u ^ E\: u{qi) = 0} has codimension 1. 
We can compute a normalized function defined to be perpendicular to this 
space Uq G {E^^. In that case we can simplify 

N 

P\{qi,x) = '^u^{qi)u^{x) = u^{qi)uQ{x). 
j=0 

If A is a 4/3-series eigenvalue born on some level mo > 1, there is an easy 
characterization of Uq; spectral decimation works "in reverse," i.e. Uq is an 
eigenfunction of Amo-i with eigenvalue i?(4/3) = 20. We can then continue 
spectral decimation back to all levels < mo because we will never encounter 
a forbidden eigenvalue. 

Theorem 3.1 (cf. [2^ Theorem 3.2, j23] Theorem 3.6). Let X be a 4/3-series 
eigenvalue with mo > 1. Then 

Am^,-iu^ = 20n^, 

i.e. Uq is an eigenfunction of Amo-i with eigenvalue 20. 

Proof. Fix a point x £ Vmo-i- First assume x is only part of a single 1-cell 
in ) and let yi,y2, ys be the other boundary points of that cell. Then 
the function v shown in Figure |5]a,is a4/3 -series eigenfunction born on level 
m (this is easy to see since the sum around any small square is 0). Now v 
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(a) * (c) 

Figure 5. Functions on ) outside of the 1-cell con- 
taining x: (a) V. (b) Coefficients in tlie linear combination 
defining w. (c) w. 



vanishes on Vmo-i \ {x}, so in particular v{qi) = 0, which forces {uq,v) = 
and hence (by Theorem 2.3) {uQ,v)mo = 0- 



Since Ug is a 4/3-series eigenfunction born on level mo, we also know that 
the sum around any small square in Tmo is 0. Taking a linear combination of 
these equations, with coefficients given by Figure [HJa, and recalling that the 
inner product weights the center 4 vertices by 2, we see that {uQ,w)mo = 
where w is given by Figure Writing {u,v + w)mo = we get 

57n^(x) + n^(yi) + u^{y2) + n^Cys) = 

or equivalently 

20n^(z) - ^o(^i) + ^o(^^) + ^o(^3) ^ 20no^(.) 
3 

i.e. Amo-iUoix) = 20uq{x). 

A similar argument works when x is instead part of two 1-cells in Fmo-i, 
with the function on in Figure [6^ playing the role of v and the one in 
Figure [6|3 playing the role of u). □ 

Another special case occurs when we fix y = qq, where qo is the center 
point of the Vicsek set. At qq, all the eigenfunctions associated with 4/3- 



series eigenvalues are equal to zero (see Section 2.2). This is a fortunate 
result because, in calculating the projection kernel at qq, all the terms from 
the 4/3-series contribute zero, so we only need to consider the eigenfunctions 
associated with 0-series eigenvalues — and these form a one- dimensional 
vector space. 



4. Numerical data for eigenvalues and eigenfunctions 

Using our implementation of spectral decimation on V5„, we can compute 
the eigenvalues of the graph Laplacians on the graph approximations 
Tm- By repeatedly applying the smallest inverse of the spectral decimation 




function R, we can use these to approximate the eigenvalues Aj of the stan- 
dard Laplacian A. Figure [7] shows plots of the eigenvalue counting function 
N(x) = : Aj < x}. Since the eigenvalue counting function grows as 
where a = log(4n — 3)/log((4n — 3)(2n — 1)), it is also useful to look at 
the Weyl ratio N{x)/x°^, shown in Figure [sj For each n, these functions are 
asymptotically periodic as a function of log x, as predicted in [T7] . What is 
rather striking and somewhat mysterious, there appears to be a convergence 
as 71 — > oo, after appropriate rescaling. We will attempt to explain some of 
this in Section [TTl 

We can also compute eigenfunctions of the graph Laplacians. Figure 
[9] shows 0-series eigenfunctions and their restrictions to the diagonal, and 
Figure 10 shows the same for some 4/3-series eigenfunctions. The eigen- 
functions in the diagonal plots have been continued with the lowest inverse 
several times to increase the number of data points. For n > 2, our 
implementation can only compute eigenfunctions restricted to the diagonal, 
./figures 11 and 12 show these plots for VS%. There is more data on the 
website ^QX- 

We observe from the data a phenomenon known as miniaturization 
Taking a 0-series eigenfunction on the mth level approximation to VS2, 
if the function is continued by spectral decimation to the {m + l)th level 
of approximation, the new eigenfunction is composed of 5 copies of the 
previous one; it is "miniaturized." Thus, eigenfunctions of higher eigenvalue 
are composed of copies of eigenfunctions of lower eigenvalue. 
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Figure 7. Eigenvalue counting function on VSn for n = 
2, ... ,7. Here m refers to the level of the graph approxima- 
tion. 



5. Spectral operators 



We can apply the theory of spectral operators to finding solutions of 
differential equations on the Vicsek set. Two major equations are the heat 
equation and the wave equation. 

5.1. Heat Kernel. The heat equation for a function u{x,t) with Neumann 
boundary conditions states 

du 

dnu{qj,t) = 0, 
u(x,0) = /(x). 



THE LAPLACIAN AND SPECTRAL OPERATORS ON THE VICSEK SET 



19 



n=2, m=7 n=3, m=5 




Figure 8. Weyl ratios on V5„ for n = 2,..., 7. Here m 
refers to the level of the graph approximation. 



Formally this is solved by 

n(x,t) = e*^/(x) 

and since the Laplacian has a discrete spectrum with an orthonormal basis 
{uj} of eigenfunctions, —Auj = XjUj, the solution to the heat equation is 

u{x,t) = 1^^"*^' (/ fiy)ujiy)My)^ uj{x). 

Usually the sum and integral can be interchanged to yield 
u{x,t) = / h{t,x,y)f{y)dfj,{y) 




6.1 728, s:0, w:31 111111111 >.:90.0081 , s:0, w:321 11111111 >.:242.591 5, s:0, w:1 31 1 1 1 1 1 1 1 1 




>.:269.3893, s:0, w:331 11111111 X:1 302.4774, s:0, w:3321 1111111 1:1 350.1 208, s:0, w:1 321 1 1 1 1 1 1 1 




Figure 9. Some 0-series eigenfunctions on VS2 and their 
restrictions to the diagonal. 
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X:2.601 8, s:1 .33, w:1 1 1 1 X:1 20.5703, s:1 .33, w:1 21 1 X:248.8754, s:1 .33, w:1 31 1 




Figure 10. Some 4/3-series eigenfunctions on VS2 and their 
restrictions to the diagonal. 
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>.:13.5167, s:0, w:31111 



>.:52.8748, s:0, w:51 1 1 1 



X:1 14.2394, s:0, w:71111 






90.3302, s:0, w:91111 



269.0855, s; 0, w: 1 1 -1 -1 -1 -1 



X: 332.3073, s: 0, w: 13-1-1-1-1 






Figure 11. Some 0-series eigenfunctions on VSg, restricted 
to the diagonal. 



X:3.1744, s:1.33, w:11111 




>.;28.1433, s:1.33, w:31111 




>.:75.7289, s:1.33, w:51111 




X:1 40.8489, s:1.33, w:71111 



>,:21 5.0306, s:1.33, w:91111 



V. 285.9038, s: 1 .33, w: 1 1 -1 -1 -1 -1 






Figure 12. Some 4/3-series eigenfunctions on VSs, re- 
stricted to the diagonal. 
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where h is defined to be 

h{t,x,y) = '^e~*^^Uj{x)uj{y) 
j 

and called the heat kernel. 

From the eigenvalues and eigenfunctions we can construct the heat kernel 
on the standard Vicsek set. This is especially easy when one of the arguments 
is the center point of the Vicsek set, since then we only need to consider 
0-series eigenfunctions. Plots of the heat kernel on the m = 4 approximating 



graph are shown in Figure 13 



Our data allows us to examine the behavior of the heat kernel h{t, qo, x) 
in greater detail. Estimates for the heat kernel are known, but they involve 
constants of unknown size. It is expected that h{t,q(),x) should involve a 
factor of t~°' multiplying a term that drops off exponentially as x moves away 



from qq. Since that data in Figure 18 suggest that the t ° factor is modified 



by an oscillating factor, we look at the ratio h{t, qo, x)/h{t, qo, qo) = H{t, x). 
Actually, it seems more plausible that 

^Jh{t,qo,qo)^/h{t,x,x) 

will be better behaved than H{t, x), but since we don't know how to compute 
h{t, x, x) effectively, this isn't an option. Note that H{t, x) is normalized so 
that H{t,qo) = 1. Also, if we ignore the influence of the boundary, which 
is certainly very slight for small t, we expect H{t/\b,Fox) should be very 



close to H(t,x). Figure 14 illustrates this invariance property. 



First we look at the behavior of H(t, x) for x restricted to the diagonal. 



Figure 15 shows some typical graphs. We also look at —\ogH{t,x), again 



shown in Figure 15 Since — \ogH{t,x) vanishes at x = qo, we try to fit a 
power law — \ogH{t,x) ~ a\x\^ where the constants a and b depend on t, 
and |x| denotes the distance to qo. However, we find that the power b varies 
significantly as we vary the neighborhood of qo where we do the fit. This 



leads us to doubt the power law model. Figure [15] shows a log-log plot of 
— log H{t, x) for some choices of t. 

There is no compelling reason to restrict x to the diagonal in studying 
the heat kernel. In a crude sense, the heat kernel h{t,qo,x) should depend 
on the distance of x to qo in the resistance metric, which coincides with 
geodesic distance in VS2- But in fact, this is not very accurate. What we 
want to look at are what might be called "heatballs" , sets of the form 

{x : h{t,qo,x) > s} 

for different choices of t and s. A naive guess would be that the heatballs 
form a 1-parameter family of sets, at least if we stay toward the center 
of VS2 where the influence of the boundary is small. Again this is only 



valid in a crude sense. Figure 16 shows some examples of heatballs for 
two different choices of t and a variety of s- values. One observation is that 
heatballs tend to spread further in directions perpendicular to the diagonal. 



24 SARAH CONSTANTIN, ROBERT S. STRICHARTZ, AND MILES WHEELER 
1=0.008333 1=0.016667 1=0.025000 




Figure 13. The heat kernel h(t,qQ,x) on VS2 for various 
values of i, on the whole Vicsek set and restricted to the 
diagonal 
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t=0.1 00000 



t=0.010000 



t=0.001000 






-0.5 -0.25 0.25 0.5 -0.5 -0.25 0.25 0.5 -0.5 -0.25 0.25 0.5 



0.6 



0.4 



0.2 



Figure 14. H{t, x) and Hit /lb, Fqx) on the diagonal (using 
the m = 4 graph approximation). 
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-log H(t,x) 



-log H(t,x) 




Figure 15. H{t, x) and — log H{t, x) for several choices of t 
(using the m = 4 graph approximation). 



Decreasing the value of s increases the size of the heatballs, so we may 
imagine that the heatballs for fixed t represent an "invasion" that spreads 
out from the center point go- By and large the invasion follows an orderly 
patter, with a cells that lies on the diagonal being invaded first at the point 
closest to go- However, there are examples where the invasion jumps around, 
and this produces examples of heatballs that are disconnected. Apparently, 
disconnected heatballs also may occur in the setting of manifolds [H] - Of 
course, it is also possible to study invasions with s fixed and t increasing. 

The trace of the heat kernel and its value at the center, when multiplied 
by f", are both periodic in \ogt (see [I3j). This is shown in ./figures 17 
and 18 on the m = 7 graph approximation. The approximate sinusoidal 



behavior is explained for the trace in [2] , and at the center in [T3] . We note 
that the approximate sines are out of phase: Fitting to a + 6sin(clogt + d) 
we get a = .90, b = .045, c = 2.33, and d = —2.2 for the trace of the heat 
kernel, and a = .4110, b = .0191, c = 2.33, and d = 1.803 for the heat kernel 
at the center. 
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t = 0.010,8 = 0.623 1 = 0.010,8 = 0.598 1 = 0.010,8 = 0.424 




5.2. Wave Propagator. The wave equation is given by 

If we impose Neumann boundary conditions and initial conditions 0) = 
at' 



0, Mzu{x, 0) = f{x), then the solution is given by 



u{x,t) = j W{t,x,y)f{y)dii{y) 
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Figure 18. Plot of f^h^t, go, Qo) versus t on the m = 7 graph 
approximation to VS2- 



where the wave propagator W{t, x, y) is given by 

W{t,x,y) = ^■'^ Uj{x)uj{y). 



From the eigenvalues and eigenfunctions we can also construct the wave 
propagator on the standard Vicsek set. As with the heat kernel, this is 
easiest to compute when on of the arguments is the center point of the 
Vicsek set, since then we only need to consider 0-series eigenfunctions. This 



is shown in Figure 19 



As already observed in the case of 5^ in [8^ the wave propagator W{t, qo, x) 
is not supported in a small neighborhood of qq for fixed t; in other words, 
waves propagate at infinite speed. This is easily explained because the dif- 
ferential operators on either side of the wave equation do not have the same 
order. However, the amount of energy that propagates at high speed is 
relatively small. So we can expect a weak substitute for finite propagation 
speed. Attempts to understand this in [8] and |6] were stymied by the com- 
plexity of the wave propagator on SG (in [2| it was shown that time integrals 
of the wave propagator are computationally tamer on SQ, but this did not 
help with a weak finite propagation speed). 

On V52 the wave propagator at the center point may be effectively com- 
puted. In particular, when we increase the level of approximation the 
graph does not change appreciably: Figure 20 shows distances between 
Wm{t, Qo, •) and Wm~i{t, Qo, •), where Wm is the level m approximation to the 
wave propagator. In Figure [19] we display the graphs for some values of t. 
Unlike the heat kernel, the wave propagator is not known to be positive, 
and indeed we see time where negative values occur. We know 



W{t,qo,x)dn{x) 



t 
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Figure 19. The wave propagator w{t, qo, x) on V1S2 for var- 
fous values of on the whole Vicsek set and restricted to the 
diagonal. 
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|w^(t,q„,.)-W3(t,q„,.)l3 



|Wg(t,q„,.)-w^(t,q(,,-)l4 



|Wg(t,qg,.)-w^(t,q„,0l5 






0.005 

0.1 0.15 0.2 0.25 0.1 0.15 0.2 0.25 0.1 0.15 0.2 0.25 



Figure 20. L?' distances between Wm{t,qo,-) and 
Wm-i{t,qo,-)- 



X unrestricted 



X restricted to the diagonal 



comparison for e = 1 e-1 2 




Figure 21. max{|3;| : \wQ{t, qo, x)\ > e} for several choices 
of e, both in the case when x is restricted to the diagonal 
and when x varies over all VS2- 



so that positive values predominate, and it seems from the data that 

\W{t,qo,x)\dfi{x) 



is bounded by a multiple of t. Recall that in Euclidean space, the singularity 
of the wave propagator worsens as the dimension increases. Our data is more 
in line with the n = 1 case. 

Our data strongly suggests an approximate finite propagation speed. We 
can quantify this by choosing a small cutoff e and looking for the maximum 
value of |x| where \wm{t, qo, x)\ > e for fixed and then letting t vary. 



In Figure 21 we show plots of this function, both in the case when x is 
restricted to the diagonal and in the case where x varies over all V52, for 
different choices of e. Notice that in both cases the slope of the function 
increases with e. 
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5.3. Spectral Projections. Another important class of spectral operators 
are the spectral projections. Let A be a subset (usually infinite) of the 
spectrum, and define 

PAf = E ^a/ 

AeA 

where 



P\f{x) =J2 ^jix)uj{y)f{y)dfi{y) 



for {uj} an orthonormal basis of the A-eigenspace. Such operators are always 
bounded on (with norm 1) and usually not bounded on or L°^. A 
natural question to ask is under what conditions is Pa bounded on for 
1 < p < oo? In the classical setting such results can be obtained from the 
Marcinkiewicz multiplier theorem [20] and we expect that analogous results 
should be valid in the fractal setting, perhaps related to the transplantation 
theorems of [TU] and [12] • We note that the results of [2T] imply that we can 
always "segment" such problems; we write A^ = A n [0, A^^] for a natural 
sequences of cutoffs that lie at the beginning of spectral gaps (in our 
case we take spectral decimation through level k). Then Pa is bounded on 
LP if and only if Paj. is uniformly (in A;) bounded on L^. 

In [2] we looked at some spectral projection on SG, but it was difficult to 
arrive at meaningful predictions because of the computational complexity 
of the data. Here we are able to examine one example in detail: the case 
that A consists of the 0-series eigenvalues. Because these eigenvalues all 
have multiplicity one, it is straightforward to compute kernels Kk{x,y) of 
the segmented projection operators Paj. for A; < 5 on VS2- We make a few 
simple observations. The first is that 



J Kk{x, y)d^{y) = 1 for every x. 



This follows from the fact that the constant 1 is in the 0-series, and every 
other 0-series eigenfunction is orthogonal to it. The second is that 

Kk{x,y) = Kj,{^{x),y) = Kk{xMv)) = Kk{Hx)My)) 

where ^ is any isometry of VS2- This is an immediate consequence of 
the fact that each 0-series eigenfunction is invariant under $ (if n is a 0- 
series eigenfunction then so is u o with the same eigenvalue, and the 
multiplicities are all one). Incidentally, we remark that invariance under all 
isometries does not characterize the 0-series spectrum; it easy to construct 
4/3-series eigenfunctions (on a higher level) that show this invariance. 

We examine the behavior of / \Kk{x,y)\dii.{y) as a function k. Table [T] 
shows the maximum over x for A; < 5. This is overwhelming evidence that 
max^. / |A'fc(x, ?/)|(i/u(y) ^ oo as A; ^ cxd, and this implies that Pa is not 
bounded in L} or L°°. Next we ask, for fixed x, what are the y values where 
\Kk[x^y)\ is large? Looking at the graphs of Kk{x,-) in Figure 22 we see 



evidence that the answer is the values of y that are close to $(x) for some 
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k 


max^ J \Kk{x,y)\dii{y) 


1 


1.4476 


2 


1.7336 


3 


2.9958 


4 


4.7955 


5 


7.6572 



Table 1. Maximum value of / \Kk{-,y)\diJ,{y) for several k. 




Figure 23. The restriction of Kk{x, •) to the diagonal when 
X is the junction point between two 1-cells, for 3 < A; < 6. 
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isometry ^. Note that for some choices of x, the set of ah <^(x) is finite, but 
for other choices it may be infinite. (For example, if x is a boundary point, 
then it is a dense subset of a Cantor set that includes the intersection of 
VS2 with the boundary of the unit square.) 

In Figure 23 we show the restriction to the diagonal of K^^x, •) when x 
is the junction point between two 1-cells, for 3 < k < 6. The behavior 
is certainly more complicated than the kernels in the standard Calderon- 
Zygmund theory. On the other hand, the graphs to appear to be converging 
to some limiting shape. It would be interesting to make this statement more 
precise, and to investigate whether there is boundedness of Pa for some 
values of p in 1 < p < 00 other than p = 2. 



6. Diagonals and the 0-series 

We can write L^(V5) = where Tio represents the eigenfunctions 

associated with the 0-series, and ^4/3 represents those associated with the 
4/3-series. These are orthogonal because the eigenvalues are distinct. 

Theorem 6.1. Each 0-series eigenfunction of the Laplacian on the Vicsek 
Set is determined by its restriction to the diagonal. 

Proof. If we look at the fractal Laplacian, we can view VS as the union of 
the diagonals with little trees attached, each tree a small copy of 1/4 V5, one 
arm of the Vicsek set. u\t satisfies —Au = \u, with dnU = at the outer 
boundary, and tt(go) is a specified value if qo is the center point, because the 
center point lies on the diagonal. 

Let v\ denote the function on 1/4VS that satisfies —Avx = Xvx, dnV\[q) = 
if (7 is a boundary point, and va(qo) = 1- To show existence and unique- 
ness, we have to show that —An = \u on 1/4 V5, 9„M(boundary) = 0, 
and u(go) = imply that u must be identically zero. Indeed, given such a 
function u, extend it by odd reflection across the center to the opposite arm 
of the Vicsek set, and set it identically zero on the other two arms. Then 
we obtain a global eigenfunction satisfying '^u{qj) = for the boundary 
points qj, so it belongs to the 4/3-series. But A is a 0-series eigenvalue, and 
by spectral decimation, there are no simultaneous 0-series and 4/3-series 
eigenvalues; the only way out is if it = 0. 

Now let T denote any tree of level m that attaches to the diagonal at y. 
Then there exists V'T : T — > 1/4 V5 with tpT{y) = Qo and '(/'T'(bdry(T)) = 
bdry(l/4V5), and 

A(/ o V'r) = (15r(A/) o ^r. 

This says that any tree can be put in one-to-one correspondence with an 
arm of the Vicsek set in such a way as to respect similarities. 
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Let u be our 0-series eigenfunction. Then (u|t) ° V't^ ~ / satisfies 
-A/ = (15)-'"A/ onl/4V5, 
dnfihdij) = 0, 

f{qo) = u{y). 

Since (15)^™A is not a 4/3-series eigenvalue (if it were, tlien so would A be) 
we have / = u{y)vi^-mx o V't- Hence 

u\t = u{y)Vl5''"^X ° V'T- 

So A and u| diagonal determine u according to the above equation. □ 

We would like to go further and say that any function in TCq is determined 
by its restriction to the diagonal, and aside from symmetry there are essen- 
tially no other conditions on the restrictions to the diagonal of TCq functions. 
We begin with the analogous statement on the discrete approximations. 

Let Dm denote the intersection of Vm with one arm of the diagonal. Note 
that = 1, = 2 and = 3^Dm-i — 1- Let Zm denote the span 

of the 0-series eigenfunctions through level m. We may consider elements 
of Zm as functions either on Vm or VS2- Note that dimZo = 1, dimZi = 2 
and dim Zm = 3(dim Zm-i) — 1, so dim Zm = #-0^ = + !)• Thus it is 
plausible that every function on Dm is the restriction of a function in Zm, 
and every function in Zm is uniquely determined by its restriction to Dm- 
In fact, these statements are equivalent. We conjecture a little more. 

Let xq,xi,X2, . . . denote the points in Dm moving from the outside inward. 
Let Dm"* = {xo, • . . , Xfc} and let Vm^ denote all points in Vm that attach to 
one of the midpoints of the intervals in Dm'^'^^ . 

Conjecture. If f ^ Zm vanishes on Dm\ then it vanishes on Vm^. In 
particular, f is determined on Vm by its values on Dm- 

This conjecture implies that there is a formula of the form 

for all / G Zm for each z E Vm and suitable coefficients cj^\ In fact f{z) = 

(z) defines the function in ^772 that takes on the values fi^Xj^ — ^jk 
Dm- The conjecture implies c~^\z) = if z S D^ . Experimental evidence 
indicates that c[.™^ can take on only 4 nonzero values, namely ±1 and ±2. 



(3) 

Some examples of are shown in Figure 24 and many more are available 



(for m = 4) on our website |7j . 

To pass from the discrete to the continuous version we consider function 
TLq n C (here C denotes the continuous functions on VS2)- Such functions 
have well-defined restrictions to D (one arm of the diagonal). To show 
that the restriction /|o of such a function determines /, it suffices to show 
that it determines f\vm foi' since Umym is dense in VS2 and / is 
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Figure 24. The values of for three choices of k. 

continuous. Let fm denote the projection of / onto Z^- By the results of 
|21] we know fm converges to / uniformly. If the conjecture is valid then 

fm'{z) = Y.c^r^ fm'{xk) for z ^ Vm and m' > m (since fm'\v„, G ^m), 

so passing to the limit f{z) = J2 ^^^k^^ fi^k) for z G Vm- Despite the fact 
that this is a finite sum for each m, it is a rather peculiar formula. The 
coefficients oscillate rapidly but do not go to zero as m increases. It does 
not seem likely that we can make any sense out of it if we do not assume 
that / is continuous. It seems unlikely that the existence of a continuous 
restriction to D for a function in J-Cq implies that it is continuous on V52. A 
more plausible conjecture is that if the restriction to D is Holder continuous 
of some order then the function is Holder continuous of the same order 
on VS2- Another reasonable conjecture is that the restrictions of TCq D C 
to D form a dense subset of the continuous functions on D. A less likely 
conjecture is that the restrictions give all continuous functions on D. 

7. Ratio Gaps 

In [5] it was shown that on SQ there exist gaps in the ratios of eigenvalues. 
As a consequence, it is possible to define operators of the form A' — a A" on 
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Figure 25. Ratio gaps for n = 2,3,4 and I = \w\ = 1,2,3. 
The gaps are the black regions. No gaps were found for 
n = 3,4 and ^ = 1. 




the product of two copies of SQ (A' and A" denote the Laplacian on each 
copy of SQ) where a lies in a gap, and these operators paradoxically behave 
in some ways like elliptic operators, despite the fact that the coefficient —a 
has the wrong sign. These operators were called quasielliptic in [5]. There 
are no analogous operators in classical PDE theory. Thus it is of great 
interest to know whether similar operators exist for products of fractals 
other than SQ. In fact [IS] shows that this is the case for VS2 and VS^. 
Also [9j investigates this question for a variant of the SQ type fractal. The 
method used in [9 , which we follow here, yields a computer-assisted proof. 
The idea is that the method introduced in [5] leads to a large number of 
tedious calculations, and these are best left to the computer. In our method 
there is a parameter £ that may be chosen at will. Increasing £ will do a 
better job finding gaps, at the cost of increasing the number of computations. 
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Let \m be a graph eigenvalue born on level mo- Then 

where Amo £ {0,4/3}. Let p be the fixed point of (/'2n-i and q = (/)2n-i(4/3). 
If Xmo = we have Am < P, while if Amo =4/3 then we either have m = tjiq 
and Am = 4/3 or Am < Q- Since q > p, we can simply write 

(7.1) Am G [0,g]U{4/3}. 

Fix ^ > 0. Any fractal eigenvalue A is of the form 

where all but finitely many of the Vj = 1. Thus there must be a word w of 
length i and some graph eigenvalue Am so that 

A = lim ° • • • ° <Ai>„„+i (Amo) 

= p™«+^+- lim pV?(</'«,(Am)) 

K— +00 

= P™»+'+>n(0»(Am)), 



^ ) . Combining this with (7.1) we see that every fractal 



eigenvalue A can be written as 

(7.2) A = /9>„(x), xE U 0^([O,g]U{4/3}). 

\w\=e 

for some integer r. 

Consider the contribution of a word ti; to the eigenvalues described by 



(7.2). If w ends in a 1, then as long as m > mo can rewrite 



o o • • • o '/'t^™o+i(Amo) = 01 o o 0*),n-i o • • • o 0t,,„(, + 1 ( Amo ) , 

for some other word w' of length £ (with one less 1 at the end) , while m = mo 
means 

4>w o 0t,,„ o • • • o (pv^^+iiXmo) = 0t«(Amo) S (0) , c/)^ (4/3) } . 



ThusjT^ is stih valid with (/>^([0, q] U {4/3}) replaced by {(py,{0),4>w{4/3)} 
in (7.2) for every word w ending in 1. Furthermore we can discard 4>w{Xmo) 
if it is forbidden. 

So far we've found finitely many intervals [oj, bi] (allowing m = bi) so that 
each eigenvalue A must satisfy 

A G p''ipn{[ai,bi]) = p''[tpn{ai),'>Pn{bi)] 

for some i and r. Therefore any ratio of eigenvalues X/ fi must satisfy 



A 



for some r, i, and j. Since pX is an eigenvalue if A is, we can restrict our 
attention to ratios A/p G [l,p] and hence to the finite number of intervals 
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p^[Rij,Sij] which intersect [l,p\. The gaps in the union of these intervals 
are then guaranteed to be ratio gaps. 



Figure 25 shows the ratio gaps that are proved to exist by this method for 
n = 2, 3, 4 using values of £ = 1, 2, 3. For all of these n there are ratio gaps 
containing ^/p^, given in Table l2j We see clearly that the number and size of 





i ratio gap 


2 3.8730 


1 [3.5370,4.2409] 

2 [3.2948,4.5526] 

3 [3.2948,4.5526] 


3 6.7082 


1 no gap 

2 [6.6952,6.7212] 

3 [6.6950,6.7214] 


4 9.5394 


1 no gap 

2 no gap 

3 [9.5357,9.5431] 



Table 2. Ratio gaps containing 
i = 1,2,3. 



fp^ for n = 2,3,4 and 



the ratio gaps increases with £. However, we have not been able to confirm 
the existence of ratio gaps for n > 5. For n = 5 none are revealed for i <2 
and our MATLAB implementation (see [7]) runs into memory problems for 
i > 3. For ^ > 3 we can, however, use a modified algorithm which searches 
only for ratio gaps containing a particular point. These searches have failed 
to find ratio gaps containing ~ 12.3693. It is not clear if these failed 
searches should be interpreted as experimental evidence for the nonexistence 
of ratio gaps, or just as evidence that we need to consider higher values of 
£ to find ratio gaps. 

8. Eigenvalue clusters 

We say the spectrum of a Laplacian exhibits spectral clustering if the 
following holds: for every integer n and e > there exists an interval / of 
length e that contains n distinct eigenvalues. 

This, for example, says you can find a million distinct eigenvalues within 
a millionth of each other. The eigenvalues will have to be very large, so it 
becomes computationally challenging to find such tight and large clusters. 
Clustering does not occur on the Sierpinski gasket SQ. Experimental evi- 
dence suggests that it does occur on the pentagasket [1] and on the Julia 
sets [TT]. The following lemma allows us to prove it holds on VS2- 

Lemma 8.1. Suppose spectral decimation holds with spectral renormaliza- 
tion factor p and spectral renormalization function R{X). Suppose R has a 
fixed point t (R{t) = t) such that \R!{t)\ > p. Then spectral clustering occurs. 
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Proof. Let (pi, (f)2, ■ ■ ■ (pN be the inverses of R{X) in increasing order. Then 
01 (0) = and 0i(O) = 1/p. There exists k such that (/>fc(t) = t and |</>'fc(t)| = 
6 < p^^, by the assumption. Choose m large enough that has at least 
n distinct eigenvalues Ai,...,A„. Then Am+j has n distinct eigenvalues 

4'k\Xi)^ • • ■ ) 4'k\Xm) and these give rise to n distinct eigenvalues 
lim p™+^'+'</.S'V?(Ap), l<p<n. 

Write g{X) = limj^oo pVi ^(A). Then 5^ is a fixed function with bounded 
derivative < M in the relevant interval of length e where we want to 

find n distinct eigenvalues. 

By taking jo large enough, we can make all the values 4'k°\^p) close 
enough to t so that < a < p'^ for all (l)-j^{Xp). This means that 

{'Pk°^''^\^p)} belongs to an interval of length no more than ca-'* where c is 
the length for ji =0. Then 

belongs to an interval of length at most cM p"^^^" (apY^ . Since ap < 1, this 
can be made < e by taking ji large enough. Thus we can find n distinct 
eigenvalues in an interval of length no more than e. □ 

On VS2, p = 15 and R{X) = 36X^ - 48X'^ + 15 A. So R{t) = t means 

2t{18t^ - 24t + 7) = with solutions 0, We are interested in the 

largest t, which is the fixed point of cp^. 

R'{X) = 3 • 36A2 - 2 • 48A + 15 = 15 + 12A(9A - 8) 

so to show R'{t) > 15 we need t > 8/9. But (4 + ^/2)/6 = 0.902... so this is 
true. We thus have clustering in VS2- Computing the largest fixed point t 
of R on VSn for n = 3, . . . , 9 we also get R'{t) > p (see Table [s]) and hence 
that spectral clustering occurs. Because the ratio R'{t)/p increases rapidly 
with n we conjecture that spectral clustering occurs for all n. 



n 


t 


P 


R'{t) 


2 


0.9024 


15 


1.6314 X 10^ 


3 


0.8905 


45 


1.3999 X 10^ 


4 


0.8891 


91 


1.2355 X 10^ 


5 


0.8889 


153 


1.1079 X 10-^ 


6 


0.8889 


231 


9.9655 X 10-^ 


7 


0.8889 


325 


8.9682 X 10^ 


8 


0.8889 


435 


8.0713 X 10*^ 


9 


0.8889 


561 


7.2641 X 10"^ 



Table 3 . The largest fixed point t of the spectral decimation 
function R on VSn satisfies R'{t) > p for n = 2, . . . , 9. 
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9. Green's Function on VSn 
The Green's function G for the Laplacian is a function satisfying 

(-AG{x,y) = 6{x,y) 

{ G{qj,y)=0 if g,- GbdryV5„. 

where S is the Dirac delta function. Then u{x) = J G{x,y)f{y)dfi{y) solves 

r -Au = /, 

\ ^Ibdry = 0. 

As a function of x, G should be harmonic in the complement of y. Suppose 
y lies in the upper right arm of V5„. The boundary points are labeled 
qi,q2,Q3,Q4, with qi corresponding to the arm where y is. Let z be the 
projection of y onto the diagonal of VSn- (In the case that y is on the 
diagonal already, z = y.) 

Now G{qj,y) = for j = 1,2,3,4. Define G{qo,y) = a, (where qo is 
the center point), G{z,y) = b, G{y,y) = c. The values a, 6, c determine G 
because G{x, y) is linear on the arms {qo, 92), (90, 93), (^o, Qa), on (go, z) and 
(z,(7i), and along the unique path joining z to y. It is constant on every 
component of the complement of these 6 sets. 

To determine the constants a, 6, c we have 3 equations that express —/S.xG{qQ, y) 
0, —/S.xG{z,y) = 0, and —AxG{y,y) = Sy. The first two equations say the 
sum of the 4 derivatives at qo (resp. z) vanish. The last says the derivative 
at y is 1. 

For simplicity assume the length of each arm is 1. (This involves rescaling 
by a factor of \/2/2 for a unit square.) Let d{z, qo) = s (so d{z, qi) = 1 — s) 
and d{x, y) = t, measured along the path. Then 

a — 6 ^ . . 
3a -\ = (at qo) 

b — a b b — c ^ , , 

+ + = (at z) 

s 1 — s t 

a — b ^ , , 
— = 1 (at y) 

Note that the third equations says c = 6 + t and the second equation says 

b — a b 
+ 



1-s 



{li z = y then t = 0, c = b, and the second equation says —AxG{y, y) 
as required.) Solving two equations for two unknowns yields 

1-s , (l-s)(3s + l) 
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In other words, 

l-d{z,qo) 
G{Qo,y) = ^ , 

^(^^y) ^ {l-d{z,qo)){3d{z,qo) + l) ^ 
G{y, y) = i^-d{z,qo)){M{z,qo) + l) ^ .^^^ 

Denote by d'{x,qo) the distance from go to the point on one of the main 
diagonals where x attaches. So d'{y,qo) = d{z,qo), for example. Then 

^ {l-d'{x,qo)){l-d'{y,qo)) 

if X and y lie on different arms. If x and y lie on the same arm and d'{x,qo) < 
d'iy,qo) then 

^ {l-d'{y,qo)){3d'{x,qo) + l) ^ 

while if d'{y,qo) < d'{x,qo) then 

Q^^ y^ ^ il-d'ix,qo))i3d'iy,qo) + l) _ 

If d{y, qo) = d{x, qo) then 

where w{x, y) is the last point on the intersection of the paths from 2; to y 
and z to x. 



10. Higher Vicsek Sets 

It is clear that VSn converges to a cross. The eigenfunctions of the Lapla- 
cian on the cross are well understood: the restriction to either diagonal is 
an eigenfunction on the unit interval, while at the center point the function 
is required to be continuous and to have the sum of its normal derivatives 
equal to zero. Thus any eigenfunction is either costt/cx on each diagonal or 
Qj sin7r(/c + l/2)x on each half diagonal, with J2 % = 0? foi' some integer A;. 
We call the first type symmetric, and the second nonsymmetric. The sym- 
metric eigenvalues (obtained by taking second derivatives) are tt^A;^, and the 
nonsymmetric eigenvalues are 7r^(A; + 1/2)^. 

We claim that the symmetric spectrum is the limit of the spectrum of the 
0-series on VSn as n — >■ 00 (these are symmetric eigenfunctions), and the 
nonsymmetric spectrum is the limit of the spectrum of the 4/3-series born 
on level 0. (the 4/3 series born on levels > 1 does not contribute to the limit 
because the eigenvalues go to infinity.) We also claim that the limits of the 
symmetric eigenfunctions are cosines, and the limits of the nonsymmetric 
eigenfunctions are sines. 
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To understand the behavior of the eigenvalues as n — > oo we can re- 
strict attention to the initial segment consisting of the 0-series eigenvalues 
Pn'>Pni'p2j-i{0)) and the 4/3-series eigenvalues p„Vn('?^2j-i(4/3)). From [21] 
we know 

r(A) = Xgn{X)K{X), and 3r(A) - 4 = /„(A)/„(A), 

so (/)fc(0) are the zeroes of gn and hn and (/)fc(4/3) are the zeroes of /„ and In- 
The zeroes of gn and fn are forbidden eigenvalues, and correspond to even 
values of k. Thus (/>3(0), (/>5(0), . . . , (/>2n-i(0) are the zeroes of hn (<^i(0) = 0, 
of course), and </)i(4/3), 03(4/3), </>2n-3(4/3) are the zeroes of In- The 
exact values 



%-i(4/3) 



2sin2f(S; 



3 3 

are computed in [23]. If j is small compared to n, which will always happen 
if we fix j and let n oo, then (p2j-ii4:/3) w ^^/6( 2n-i )^' Since tl^n{t) ~ t 
for t near we have 

/3nV'n(</'2,-i(4/3)) ~ (2j - 1)' - ^(i - 1/2)' 

as n — > 00. 

There is no exact computation of the zeroes of hn, but the zeroes of gn 
are known, so 

l + cos(^^)vr 2sin2|(^^ 



fn\ ■ ^ 2n-l /" ^'^^^^ 2V2n-l. 
2j(Uj = 



3 3 
and we have interlacing of zeroes of gn and so 

</'2i-2(0) < 02j-l(O) < 02i(O). 

This implies 

4n — 3 TT^ , ,9 , ,, , 4n — 3 vr^ , ,9 

—-(2, - 2f < pnM^2,^m) < (2^)^- 

If we assume that the lower bound is the asymptotically correct value, then 
we obtain the expected value ^^(j — 1)' for the limit. We will show below 
that this is indeed correct. 

We can also understand why the VSn eigenfunctions, restricted to the 
cross, converge to the eigenfunctions of the cross. To see this, we look at the 
graph eigenvalue equation on Vi . Note that Vi consists of four arms of n — 1 
squares joined at a central square. We label the diagonal vertices of one 
arm xi, X2, . . . , x„ and the below and above diagonal vertices yi, . . . , i/n-i 



and zi, . . . , Zn-i (see Figure 26). By symmetry we will have u{yj) = u{zj) 
for every eigenfunction. The eigenvalue equation (with eigenvalue Ai) at 
says 

(1 - \i)u{yj) = \{u{xj) + u{xj+i) + u{yj)) 
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Figure 26. Vertices on one arm of Vi. 

So we obtain 

u{Xj) + u{Xj+l) 

u{zj) = u{yj) = 2-3Ai 
For 2 < j < n — 1 the eigenvalue equation at Xj is 

(1 - Xi)u{xj) = ^(u{xj-i) + u{xj+i) + 2u{yj) + 2u{yj+i)^ 

If N , / N , 2u{xj-i) + 4u{xj) + 2u{xj+i)\ 

We can simplify this equation to 

(1 - 3Ai)7i(.T,-) = i(n(xj_i) + u{xj+i)). 

Note that this is exactly the eigenvalue equation (for eigenvalue 3Ai ) on the 
interior of the linear graph xi,...,Xn- Similarly, at the endpoint Xn the 
eigenvalue equation is 

(1 - Ai)n(x„) = ^(^n(x„_i) + 2u{yn-ij^ 

1 / 2u{Xn-l) + 2u{Xn)\ 

= 3 r^'^--'^ + 2^^ j ' 

which simplifies to 

(1 - 3Xi)u{Xn) = u{Xn-l), 

and this is the correct eigenvalue equation (for eigenvalue 3Ai) with Neu- 
mann conditions at that endpoint. 

The equation at the endpoint xi will depend on whether we are looking 
at the 0-series or the 4/3-series. For the 0-series the values along all four 
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arms will be identical, so the eigenvalue equation is 



(1 - Xi){u{xi)) = ^(3n(xi) + u{x2) + 2n(yi) 



' 3n(xi) + u{x2) + 



6 V 2 - 3Ai 

which simplifies to 



{l-18^-^^X?ju{xi) = u{x2). 



For the 4/3-series the sum of the values on all four arms will be zero, so the 
eigenvalue equation is 



Ai)n(xi) = ^( - u{xi) + u{x2) + 2u{yi)'^ 



6 



6 V ' ' ' ' 2 -3Ai 



which simplifies to 



(3 — 6Ai)m(3;i) = u(x2). 

These should be compared with the eigenvalue equation for the eigenfunc- 
tion u with eigenvalue 3Ai on two copies of the linear graph with even and 
odd symmetries, namely 

(1 - 3Ai)'u(xi) = l{u{x2) ± u{xi)). 

Note that we get the identical equation in the odd case, but in the even case 
we get 

(1 - 6Ai)n(2;i) = u{x2), 

so there is a significant distinction. In the case of the 4/3-series, we can 
therefore identify the restriction of the eigenfunctions to the diagonal with 

/2k — 1 \ 

u{xk) = sin7r(i - 1/2) ( 2n - 1 ) ' 1 < J < »^ - 1- 



Figure 27 shows some 0-series eigenfunctions plotted against the symmet- 
ric eigenfunctions on the cross for n = 3, 6, 9. It appears that u closely 
approximates 



n(x.) = cosvrj(?^) 



We now sketch a proof that the eigenvalues 3Ai approach 



3Ai = 1 — cos = 2 sin 



2n - 1 2n - 1 

and the eigenvectors u{xi^) approach u{xk) as n ^ oo. Here we fix the value 
of J, and we require the appropriate error estimate since both 3Ai and 3Ai 
tend to zero. The idea is to use standard perturbation theory, using the fact 
that the two eigenvalue equations differ only at the single point xi, and the 
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n=3, >.:13.715, w:3111 



n=6, >.:13.606, w:3111 



n=9, >i:13.4831, w:3111 
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-0.5 
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n=3, >.;53.728, w:5111 



n=6, X;52.1354, w:5111 



n=9, >.;53.0116, w:5111 






n=3, J.:279.1452,w:5211 



n=6, 07.7091, w:71 11 



n=9, >.:1 15.6708, w:7111 






Figure 27. 0-series eigenfunctions on V<S„ for n = 3,6,9 
plotted against the corresponding symmetric eigenfunctions 
on the cross. 



fact that the eigenvector u{xk) is fairly uniformly distributed, so the value 

u{xi) is relatively small. 

Let E denote the symmetric n x n matrix of tridiagonal form, with 



Ekk 



1 2<k<n-l 
^ /c = 1 or n 



and E^i^^i-^ = E^]^_i-^ = —1/2. Let G denote the diagonal matrix with 



Gkk 



1 1 < A; < n- 1 
h k = n 
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and let G denote the diagonal matrix with 

{1 2<k<n-l 
\ k = n 

3Ai-3 u _ 1 
3Ai-4 «^ — 

Then the two eigenvalue equations may be written 

Eu = 3AiGn, 
Eu = 3AiGn. 

Note that the first equation is not a linear generalized eigenvalue equation 
because G depends on Ai, but this does not really matter in our argument. 

The gist of the argument is that G — G '\s a matrix with only one non-zero 
entry (Gu — Gu) and we can bound this entry since Ai is bounded away 
from 4/3 for the 0-series; and also we know u exactly, hence |u(xi)| < 1 
while {Gu,u) = n/2. This yields the estimate 

{{G-G)u,u) 

(Gn,n) U 

With a little more work, we can get the estimate 

{{G-G)u,u) ^^(l 



{Gu, u) 



n 



for the first N eigenfunctions (A^ is fixed as n — > oo). But this is exactly 
what we need to estimate Ai — Ai. If we take the inner product of the first 
eigenvalue equation with u and the second eigenvalue equation with u, then 
using the symmetry of all the matrices we obtain 

\i{Gu,u) = Xi{Gu,u) 

hence 

A.-A, = A,i(^^^ = 0f4 
{Gu,u) W 

for the first N eigenvalues, since we know Ai = 0(l/n2). With a little more 
work we can show that u — u = 0(l/n) when u is properly normalized. 

So far we have dealt with the level 1 eigenvalues Ai. The actual eigenvalues 
A on VSn are given by A = ipniM) for the lowest segment of the spectrum 



(this will include the first eigenvalues once n is large enough). Figure 28 
gives experimental evidence for the estimate t < ipnit) < t + ct^ on < t < 1 
for a constant c independent of n. This shows that A — Ai = 0(l/n^) as 
n — > oo for the first N eigenvalues. 
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v|( (x)/x for n=2,...,12 




Figure 28. tl}n{x)/x on [0, 1] for n = 2, . . . , 12 (V'n+i > V'n)- 
11. Weyl Ratio 

We now describe in more detail the Weyl ratio Wn{t) = Nn{t)/t°'"' on 
VSn, where 

log(4n — 3) 



log Pn 



and 



Nn{t) 



E 



m{Xj 



is the counting function for the number (counting multiplicity) of eigenval- 
ues. According to a general theorem of Kigami and Lapidus |17!, Wn{t) = 
limfc_»oo Wn{Pnt) exists. In order to compare Wn for different values of n, we 
normalize by Wn{s) = tt;„(Aip^) so that Wn is a periodic function of period 
1 with tt)n(O) = t(;„(Ai). 

From the data it appears that Wn is converging to a limit as n — > oo, 
but this limit has nothing to do with the Weyl ratio on the cross, which 
tends to a constant. While we cannot supply a complete explanation of this 
phenomenon, we can make a few observations about the behavior of Wn{s) 
for some values of s. Because of high multiplicities the functions Wn and Wk 
have jump discontinuities. We write Wn{t~) = lim^^^- Wn{s) and similarly 
for Wk{s~). First we note that it is possible to compute Wn{\i) for small 
values of i. 



Lemma 11.1. For i < j < n — 1 we have 

^ (A2i-l)"" 



Wn{X2j~l) 



1 



(A2i-l)"" 

Proof. A simple induction argument shows that 

iV(pJ^A2,-i) = (4j-l)(4n-3)'= + l. 



4j 



(A2,)°^ 



m(p„A2j- 



2(4n-3)'' + l, 
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since p^X2j-i is a 4/3-scrics eigenvalue born on level k. Thus 

, ^ (4i - l)(4n - 3)fe + 1 ^ (4j - l)(4n - 3)^' + 1 

'^"'^ (ft';A2,_i)"" (4n - 3)^(A2,_i)"'' 

and the computation of Wn{Mj-i) follows by taking the limit. Similarly 
we obtain the result for 'u;„(A^_]^). We can also show by induction that 
N{p^X2j) = Aj{An - 3)*^ + 1, and the result for WniMj) follows. □ 

In particular, we have Wn{^~) = 1/A°" and Wn{^) = 3/ A"". As n ^ 
oo we have Ai — > 7r^/3 and Oi^ — > 1/2 so liniy^ — >oo 

Wn{0 ) = \/3/tt and 

lim^^oo ^i'n(O) = sVs/tt. It is more difficult to get information about lim- 
iting behavior of Wnis) for other values of s because we would have to 
simultaneously let j — cx) as n — oo. Although we know A; —>■ |(^)^ for 
fixed / as n ^ oo, the convergence is not uniform in /. 

Because the spectrum has large gaps on either side of p^Ai, we can 
say more about the behavior of Wn{s) for s near 0. In fact the eigen- 
value just below p^Xi is PnV'n('?!'2n-i(0))> the eigenvalue just above it is 
p^Vn(.^241i(0)). So for 

PnM^S-im < t < p^Ai, 

the value of Nn{t) is (4n - 3)'' so Wn{t) = (4n - 3)''7t"". In taking the 
limit as A; — ^ oo we note that ^2n-i(0) ~^ Pn-, the fixed point of (j)2n-ii so 
Wn{t) = for ipn{Pn) <t < \\ ov equivalent Iv 

~/x 1 f /log Ai - log Vn(Pn)\ , 
^n(«) = for - — < S < 0. 



Ar(/>^)- V logPn 

Similarly Wn{t) = 3/i"" for Ai < i < pnfpni4'2{Pn)) , or equivalcntly 

3 f f log pn + log IpniMPn)) -^Og\l\ 

'Wnis) = / ,v N for < s < ^ — — . 

^ ^ A°"(p^")-- - - V log Pn J 
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